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Abstract 

A 6 term, which couples to topological charge, is added to the lattice 
CP^~^ model. The strong- coupling character expansion is developed. 
The series for the free energy and mass gap are respectively computed 
to tenth order and fourth order. Several features of the strong-coupling 
analysis emerge. One is the loss of super confinement. Another is that 
in the intermediate coupling constant region, there are indications of 
a transition to a deconfining phase when 9 is sufficiently large. The 
transition is like the one which has been observed in Monte Carlo sim- 
ulations of a similar lattice CP^~^ action. 
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I. Introduction 

One of the unresolved puzzles of the standard model is the strong CP problem. 
The problem emerged after the discovery of instantons in four-dimensional Yang- 
Mills theories Instantons represent barrier penetration processes between different 
classical n- vacua. The "true" quantum vacua are believed to be linear combinations of 
n- vacua weighted by a phase exp (inO) [0, |^. Alternatively, one can add 6 J df'xv{x), 
where v{x) = g'^F°'F°'{x) / (327r^), to the action lagrangian. Here, / d'^xv{x) = Q is the 
topological charge. It is one (minus one) for an instanton (anti-instanton). The action 
term 9 J d'^xv{x) breaks parity, time-reversal invariance and CP symmetry. Although 
the ^-term is a total derivative so that it does not contribute in perturbative Feynman 
diagrams, instantons and topological quantum fluctuations contribute to it. Such 
topological fluctuations are known to exist because they contribute significantly to 
the 7]' mass Lattice simulations have computed these topological fluctuations 

and yield results of the right magnitude to explain m^/ P] 



When quarks are included, phases in the quark mass matrix M. contribute to CP 
violation also. Using a U{1) axial rotation, a common phase in M. can be shifted to 



the 6'-term due to the axial anomaly [|14| , [151 . The physical effective theta parameter 
9eff is 6'e// = 9 + ArgDetM.. A non-zero value of 6'e// contributes, among things, to 
the electric dipole moment of the neutron. Requiring compatibility with experiment 
leads to 9eff < 10^^ The difficulty in understanding how can naturally 

be so small constitutes the strong CP problem. 



Various solutions to the strong CP problem have been proposed |jT8|, |T9|, ^ 



Perhaps the most theoretically attractive suggestion is the Peccei-Quinn mechanism 



It predicts the existence of an axion [p2| , p3| . Such an axion has not been seen in 
accelerator experiments. A modification of the mechanism leads to an invisible axion 
Astrophysical considerations and accelerator experiments have left open a 
somewhat narrow mass window for the invisible axion pH|. 

Another possible solution is that the current- algebra mass of the up quark is zero 
3^, Various analyses suggest that the up-quark mass is non-zero and equal 



to about 5 MeV [^, |3^. However, the up-quark mass has been the subject of much 
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debate and some theorists feel that a zero value is not ruled out ISi 



The Peccei-Quinn mechanism involves a relaxation field which allows 6'e// to be 



rotated away. This field is the axion. In ref. [Q, a 2 + 1 dimensional model was 
considered in which the Yang-Mills sector generated its own relaxation field. This led 
to a natural solution to the strong CP problem within the pure Yang-Mills sector. 
Criteria were established to determine when such a natural relaxation mechanism 
arises [^. Unfortunately, it has not been possible to determine whether or not these 
criteria are satisfied for four-dimensional Yang-Mills theories. 

The strong CP problem in a pure Yang-Mills theory involves vacuum physics. 
Vacuum physics is governed by the long-distance behavior of a theory. The long- 
distance behavior of non-abelian gauge theories involves non-perturbative and strong- 
coupling effects. Since this regime is not so theoretically well understood, there is the 
possibility that some subtle mechanism exists which resolves the strong CP problem 
naturally. 

Given the complexity of the strong CP problem, it is useful to consider simpler 
systems. One class of such systems is the two-dimensional CP^~^ models. These 
models are asymptotically free and have instantons for all > 2. In addition, 
a ^-term can be added to the action. The analysis in the large limit leads to 
considerable insight into their physics [^, |^. The 6 dependence is expected to be a 
effect - this contrasts the exponentially suppressed behavior exp(— cA^) expected 
from an instanton analysis ||^. In the leading order result, the CP^^^ models 
consist of free particles and A^ free antiparticles, and there is no 9 dependence. 
When the first 1/A^ correction is taken into account, quantum fiuctuations produce 
a dynamically generated U{1) gauge field which, in turn, generates a constant force 
between particles and antiparticles, thereby leading to confinement. The free energy 

per unit volume behaves as cO'^/N, where c is a constant. It is expected that 
the confining force becomes stronger and stronger as higher-order 1/A^ corrections 
are included since the CP^^^ models exhibit super confinement [^]. In a theory 
with superconfinement, the confining force is so strong that an antiparticle must be 
positioned on a particle. Not only is there confinement at asymptotic distances but 
there is local confinement as well - particles and antiparticles in a bound state cannot 
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be separated at all. The dramatic changes in behavior in going from leading order 
to next-to- leading order are indicative of the singular nature of the expansion 
in the CP^^^ models. The singular nature is also reflected in the fact that strong 
coupling and expansions do not commute |^ . 



G. Schierholz and co-workers have performed Monte Carlo simulations on the two- 
dimensional CP^ model with a 9 term [^. A priori surprising results have emerged. 
For fixed inverse coupling (3, a phase transition occurs as a function of 9. Letting 9(. 
denote the critical value, the free energy per unit volume is well represented by 

^_{a{[5)9^ 9<9, 

\c{^) 9>9, ■ ^'-'^ 

In other words, for 9 > 9^, the free energy has no dependence on 9. In two dimensions, 
this implies that the string tension vanishes for 9 > 9c so that confinement is lost. 
Furthermore, 9(. appears to go to zero as the coupling g goes to zero (or (3 goes to 
infinity). The behavior of for ^ < 6'c is compatible with the next-to- leading order 
result. This is not the case for 9 > 9c. Since additional corrections have 
not been computed, it is unclear whether there is a discrepancy. 

The two-dimensional Schwinger model has a cusp in the free energy at 9 = n, 
signaling the spontaneous breaking of CP invariance [Q. A similar phenomenon 
is expected in the CP^^^ models. In fact, at infinite coupling, j3 = 0, the lattice 
(jpN-i models do undergo spontaneous CP breaking at = vr Hence, a phase 
transition with = is not unexpected. The surprising feature of the CP^ model 
is the dependence of 9c on (3 It suggests that to obtain a confining theory from 
the continuum limit of a lattice CP^^^ model, one must take 9 to zero. Thus it 
appears that continuum confining CP^~^ models must have 9 = Q. If the analog 
statement were true for a four- dimensional Yang-Mills theory, then the strong CP 
problem would be solved. 

An analysis of the free energy of the four- dimensional SU{2) Yang-Mills the- 
ory indicates that the free energy behaves as in equation (|1.1| ). In contrast to two- 
dimensional systems, however, a constant JF behavior for 9 > 9c does not necessarily 
indicate a loss of confinement. Nonetheless, if 9 must be selected to be less than 9c for 
certain as-of-yet-unknown physical reasons and 9c goes to zero as P ^ oo, then the 
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strong CP problem would be solved. Up to now, simulations of the four-dimensional 
SU{2) gauge theory have been done at only one value of /3, so that it is not known 
whether 9c goes to zero as /3 ^ oo. 

Since the strong CP problem is related to long-distance physics and since long 
distances are equivalent to strong couplings in an asymptotically free theory, it is 
of interest to perform lattice strong-coupling expansions. As far as we know, only 
one strong-coupling expansion has been performed in the presence of a 6' term. The 
zeroth order and order terms for the free energy and mass gap were obtained in 
ref. [^]. Strong-coupling series are best organized in character expansions. Char- 



acter expansions were introduced in ref. for lattice gauge theories. They were 



obtained for matrix- model spin systems and for the large N limit in [Q]. Finally, the 
adaptation of such methods to vector-like spin systems, including the CP^^^ models. 



was achieved in [^, ^ . As far as we know, the techniques for performing character 
expansions in the presence of a ^ term have not yet been developed. One of the 
purposes of the current work is to fill this need. Another purpose is to gain insight 
from such analytic strong-coupling computations into the physics associated with 9. 
Strong-coupling analysis is often one of the best ways to gain an understanding into 
the non-perturbative long-distance behavior of asymptotically free theories. It was 
by these means, for example, that K. Wilson demonstrated that non-abelian gauge 



theories have the potential to confine quarks in four dimensions ^ . 

II. Methodology 

The action for the continuum CP^~^ model is 

S = (3N J d'^x (d^z*d^z' + (z*d^z^ ' ^^-^^ 

where z{x) is an A^-component complex field satisfying z*z^{x) = 1. Eq. pos- 
sesses the local symmetry z{x) —>■ exp{ia{x))z{x). Hence, the local gauge-invariant 
degrees of freedom at a point x live on the manifold CP^~^. Using an auxiliary field 
A^{x), the action in Eq. (|2.1|) can be written as 



S = (3N d^x {d^ - iA^,) z* {d^ + z' . (2.2) 
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Quantum fluctuations turn A^{x) into a propagating field. It becomes the U{1) field 
mentioned in the Introduction and confines the charged +1 particles z* to the charged 
— 1 particles Zj. 

There are always many equally satisfactory ways to latticize a continuum action. 
For the CP^^^ models, two lattice actions naturally arise - they correspond to the 
naive discretization of Eqs. ( [^.ID and ( |2.2| ). A lattice version of ( |2.1| ) is 

S = f3NY^ ^■+A ■ Zx zl- z^+A , (2.3) 
while a lattice version of ( p.2|) is 

S = f3N ^ i^zl ■ z^+A exp (^v4^,x+a) + ■ exp (-iA^,x+A)) • (2.4) 

a;,A 

In Eqs. ( ^.31 ) and ( |2.4| ), the sum over A involves the c? positively-directed nearest 
neighbors to x, so that A takes on the values ei, 62,. . ., e^, where e, is a unit vector 
in the ith direction. In what follows, we restrict to the case of two dimensions so that 
d = 2. 



For the action in ( p.3|) , a 6 term is defined as follows |^8|. Associate a 11(1) field 
[/(x, x + A) with each link via U {x,x + /S) = z*{x) ■ Zx+a/\Zx{x) ■ z^+aI- Let Up 
be the product of the U {x,x + A) around a plaquette. Define the local topological 
density Up via Up = log (Up) / (27r). The branch of the logarithm is taken to be between 
— vr and tt. The total topological charge Q is Q = J2p ^p- The 6 term to be added to 
the action is iOQ, that is, 

i9 

Se term = —J^ ^'^9 (Up) . (2.5) 
271 p 

The 6 term is a complicated function of the z{x). 

For the auxiliary gauge field formulation in Eq. (|2.2| ), a term can be defined by 
setting the link variable U {x,x + A) to be exp (iA^^x+A), letting Up be the product 
of the link variables around a plaquette and using Eq. (|2.5| ). 

At strong coupling, field configurations fluctuate considerably and might lead to 
bad behavior in Eq. (|2.5|). However, fluctuations are not so severe as to render Eq. 



( p.5|) so singular that strong-coupling computations cannot be performed. In fact. 
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ref. found no computational difficulties in leading order. Likewise, no difficulties 
arose in the higher order calculations of the current work. 

Boundary conditions determine whether Q is quantized. The analysis can be 
performed in a model-independent way. Since one is dealing with gauge-invariant 
functions of a compact U{1) field, one can transform to a set of independent gauge- 
invariant degrees of freedom. In two- dimensions with open boundary conditions, the 
plaquette variables Up are such a set. Writing Up = exp{ifij) with — tt < fij < n, 
Figure 1 shows a configuration of U{1) fields which produces the most general con- 
figuration of the Up. In this figure, most link variables equal 1. A Wilson loop which 
runs around the boundary of the lattice has the value exp(i fij)- Since this factor 
is not necessary 1, Q need not be an integer. 

With periodic boundary conditions, the most general configuration is that of Fig- 
ure 1 with exp{iJ2ij fij) = 1 and for which all horizontal links are multiplied by 
exp(iai) and all vertical links are multiplied by exp{ia2). These constant phases 
do not change the value of Up but allow the Polyakov lines at the lower-most and 
left-most sides of Figure 1 to be non-zero. Since the condition exp{iJ2ij fij) = 1 
implies that fij = 27iQ, Q is integer. The constraint can be handled by inserting 
J2q S {j2ij /ii/(27r) — into the functional integral. This makes the dependence on 
6 periodic since the following combination arises in the functional integral 

exp [i^ E E S /.,/(27r) - Q j = E exp [i^-^^ E /^.) • (2-6) 

Here, the resummation formula J2q ^ (j2ij /ii/(27r) — = Sm^xp (imj^ij fij) has 
been used. 

For the action in Eq. ( p. 41) , the change from link variables to plaquette variables 
Up can be explicitly performed with unit Jacobian. At /3 = 0, the integrand for the 
partition function Z{[3 = 0) is Eq. (p.6|) , where the measure is Yiij f-j^dfij/ {21:) 

'0 + 27rm' 



Z(/3)|^=o = E 



sm 



V 



(2.7) 



e + 27rm 

where V , the number of sites of the lattice, is the volume of the system. The free 
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energy T per unit volume, defined by JF = — limy^oo(log^)/^, is 

T = -\o^{^^€v^{^ , for-7r<^<7r , (2.8) 

and T is determined by periodicity in Q when Q is outside the range from — tt to vr. 
Eq. (|2.8| ) is the result obtained in ref . |^| . In what follows, we use periodic boundary 
conditions so that all quantities are periodic functions of Q with period 27r. Formulas 
below are written for the case —7i<6<tt. 

To develop strong-coupling character expansions, we have found it calculationally 
more convenient to use the auxiliary gauge field formulation. Hence, throughout 
the rest of this work, we use the action which is the sum of Eqs. ( |2.4| ) and ( |2.5| ). 
The expansion procedure is similar to the one discussed in ref. One performs 

a group-theoretic-like Fourier decomposition of the action associated with each link. 
However, because of the 6 term, one cannot immediately perform the integrations 
over the auxiliary gauge fields. One writes 

exp (3N [2;* ■ z^+A exp {iA^,x+A) + ■ ^x+a exp (-^A^^^+a) = 

^o(/3) d(^m;n) exp [i{n - m)v4^.,^+A] 2;(m;„)(/5) /( m;n) l^^x; ^x+A ) , (2.9) 
m,n 

where Zq{(3) is a normalization factor and where the sum is over non- negative integers 
m and n. Character-like "representations" r are governed by the pair {m\n)\ r = 
{m;n). Here, dr, Zr and fr are respectively character-like dimensions, expansion 
coefficients and representations for the U{N) vector models See Appendix A for 
explicit formulas for low orders. 

With the appropriate normalization of the f\r) , one has 



/ 



dv fr{w, V) f*,{w', ^) = ^ ^rr' fr{w, w') , (2.10) 



and 

fr{w, w) = dr . (2-11) 

The calculational procedure is to expand the action of each link as in Eq. ( |2.9| ) and 
then perform the integrations over the z{x) using Eq. ( p.lOj ). Lastly, one does the 
integration over the auxiliary gauge fields. Such integrations involve Wilson loops. 
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/ dUpexp 


'ilog {Up) 


U " 
t^p 


J dUpexp 


'^log {Up) 





Using Stoke's theorem, such Wilson loops can be written as products over plaquettes 
of plaquette variables to various powers. The result 



(-1)" (2.12) 

is useful. 
III. Results 

It is straightforward to compute the partition function using the character-like 
expansion in Eq. ( |2.9| ) and the above-mentioned procedure. Exponentiation of volume 
factors naturally occurs so that Z = exp(— V^JF) where is the free energy per unit 
volume. We have computed to tenth order, which is four more terms than in ref. 
|1[. When expanded in f3 up to order our results agree with those of ref. [El 



The fourth, sixth, eighth and tenth orders provide new information about the CP^~^ 
models in the presence of a 6' term. Because the result for JF is rather lengthy, we 
have relegated it to Appendix B. See Eq. (p.l| ). 

An analysis of the strong-coupling expansion leads to some immediate physical 
results. The first is the loss of super confinement when 6^0. When 6^0, strong- 
coupling diagrams appear in fourth and higher order terms in which a positively 
charged loop does not have to be directly on top of a negatively charged loop. Local 
charge liberation occurs. Charges are no longer super confined. Instead they are 
bound together by ordinary confinement. Thus, the confining force is weakened. 
Furthermore, charged loops oriented counterclockwise arise with probabilities which 
are different from loops oriented clockwise. This is expected and is an indication of 
the explicit CP violation generated by a term when 6 ^ or 6 ^ n. 

For (3 small, higher order terms do not significantly change the result in Eq. ( |2.8| ): 
JF increases as |^| increases and its graph is concave upward. As f3 becomes larger, JF 
is lowered particularly for larger values of 6. Then, for sufficiently large (3 a maximum 
in JF arises at 6 near ivr. Figure 2a illustrates this for the = 4 case. Other values 
of A^ are qualitatively similar. As j3 is further increased, the maximum in JF occurs 
at smaller values of 6. See Figure 2b which plots JF versus 6 when (3 = 0.7. For the 
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= 4 case, the maximum eventually reaches 9 = for /3 ^ 0.90. However, at these 
moderate values of (3, higher order terms might be important. 

In two dimensions the string tension a{e,6,l3) between a charge e and a charge 
— e can be obtained from the free energy via the standard formula cr(e, 6, (3) = J-'{6 + 
27ie,P) — Hence, the value of 6 at which the free energy stops increasing 

and begins to decrease marks the point where confinement is lost for infinitesimal 
charges. Using our results for JF, a phase diagram can be constructed. It is shown in 
Figure 3 for the CP^ model. The reader should be forewarned that this phase diagram 
is obtained under the assumption that strong coupling results can be extrapolated 
to moderate values of /3. The dotted line in Figure 3 indicates where CP invariance 
spontaneously breaks. This line bifurcates at around j3 = 0.56. For (3 larger than 0.56 
there are two phases as a function of 6. The region labeled "infinitesimal confinement 
loss" is where charges of very small magnitude no longer experience a confining force. 
The phase diagram of Figure 3, which is obtained by analytic methods, is qualitatively 
similar to the diagram of refs. [^, obtained from Monte Carlo simulations of the 
adjoint form of the CP^ action in Eq. ( |2.3| ). We have also taken our tenth order 
series, expanded it in j3, and obtained Fade approximates to try to increase the range 
of validity of our strong-coupling results. For the diagonal approximate, the f3 at 
which infinitesimal confinement loss first takes place is shifted to larger values. Even 
for f] as large as 2.0, there is confinement for 6 small. Hence, the phase diagram using 
the diagonal Fade approximate looks as in Figure 3 but with the lines separating the 
regions "full confinement" from "infinitesimal confinement loss" shifted to the right. 

The mass gap m can be extracted from the two-point function: 

J2{z{0, 0) ■ z*{x, L) 2*(0, 0) ■ z{x, L)) ~ exp(-mL) (3.1) 

X 

Summing over x to project onto zero momentum makes it easier to compute m. The 
strong-coupling series for m at 6* = has been obtained to tenth order in ref. . 
The computation for 7^ is considerably more complicated in that sets of graphs 
must be summed. Here, we compute m to order four when a 6 term is present. This 



is one more order than the j3 expansion computation in ref. [41 1. Our results agree 



with [it] and [41 1 for those terms which are in common. The computation of m is 
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lengthy and is relegated to Appendix B. See Eq. ( |B.2| ). 

The second order ^-dependent term lowers the value of m. This is physically 
expected. The loss of superconfinement means that the charges can move more freely 
within the bound state, thereby lowering their constituent mass. However, the fourth 
order ^-dependent term counteracts the effect and tends to increase the value of m 
particularly for large 6. In Figure 4, we plot m versus (3 for various fixed values of 
6. For sufficiently small (3, m is lowered when 6 is non-zero, but the effect is small 
and not too visible for the scales used in Figure 4. The curves for < 1.0 are almost 
identical to the case ^ = and are displayed as one curve. For 6 > 1.0, fourth order 
contributions increase m significantly at moderate values of /?, however, at such (3 
values, higher order terms can be important. 



For sufficient large /?, refs. |^ has observed a phase transition as 9 is increased. 
A priori such a transition could be due to the vanishing of m. Our computation of m 
suggests that this is not the case and supports the conclusion that the transition is 
due to the loss of confinement. If confinement is lost, one expects on physical grounds 
that m = 2mo > where ttiq is the mass of a free charged particle. 
IV. Conclusions 

In this work, we have developed and performed the first character expansions for a 
lattice model with a 6 term. The purpose was to gain insight into the strong-coupling 
behavior of the physics associated with a topological term. 

We found the following physical effects. The superconfinement property of the 
(jpN~i jg jQg^ £qj^ 0^0. When ^ 7^ 0, superconfinement liberation occurs. A bound 
state still persists but is tied together by ordinary confining forces. The weakening of 
the binding in the bound state is refiected by the lowering of its mass at small (3. 

In Monte Carlo simulations, refs. observed a phase transition when 9 was 
increased. It was argued that this transition was associated with the loss of con- 
finement. Our analytic computations, although not conclusive and performed with a 
different lattice version of the CP^~^ model, support this conclusion. Furthermore, 
the mass of the bound state does not appear to vanish at the deconfinement point. 

The main motivation for our work is to gain insight into the strong CP problem. 
Our analytic computations support the conclusion of refs. |^ that the analog of 



10 



the strong CP problem is naturally resolved in the CP^~^ models. It was argued 
that, as /5 — i> 0, 6' must be adjusted to as the continuum limit is taken in order 
to have a confining theory with bound states. Our analytic computations support 
this statement. Assuming that higher-order corrections do not invalidate results, the 
strong coupling computations show that 9 must be decreased as /? is increased to 
remain in the confining phase, but they are unable to say whether 6 must be taken 
to zero as /? — oo. 

The conclusion concerning the strong CP problem appears rather specific to two 
dimensions. In two dimensions in the continuum, a 6 term in the action produces 
an immediately physical effect - it introduces charges at spatial infinity which create 
a constant background electric field. This has a direct effect on a bound state. In 
four- dimensional gauge theories, a 6 term affects physics more indirectly: If a theory 
possesses monopoles, then such monopoles acquire an electric charge proportional to 
6 ||50|. In Monte Carlo simulations at a particular value of the coupling, ref. |^ 



has shown that the SU{2) four- dimensional Yang-Mills theory has a dependence on 
6 similar to that of the CP^~^ models. In particular, a phase transition occurs at a 
certain value of 9. Given the difference between two and four dimensions, it is not 
probable that this transition is related to deconfinement. It is thus quite important 
to determine the physics associated with the transition to decide whether 9 must be 
tuned to zero as the continuum limit is taken. 
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Appendix A: Characterlike Expansion Quantities 

The leading order character-like representations used in our computation are 



l N{N+l) ^ , 

f{2-fi){.V,w) = \ [V -W) , 
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f , w 



lN+1 



fi2;l){v,w) 
f(m;l)iv,w) 



, N-1 

{N+1 



V ■ wv ■ w 



1 

N 



(A.l) 



V ■ wiv ■ w ) — 



\ 2(iV-l) 



N+1 



V ■ w 



and the corresponding dimensions dr are 



i(0;0) 



d{i-o) = VN , 



d 



{2;0) 



'N{N + 1) 



d(i 



(i;i) 



dil: 



{l;m) 



d(m;l) ■ 



^N{N -1){N + 2) 



(A.2) 



These formulas can be derived with the help of the integral 

(iV-1)! 



dvvi^ Vi^...Vi^v*^v*^...v*^ 



V 6'^ 5'^ ...6' 



(A.3) 



where the sum is to be taken over all permutations of the j indices. In Eq. (|A.3|) and 



throughout this appendix we are using the following normalized notation for vector 
integrals 

/*,,...). . (A.4) 

7 \dvb{v*-v-V) 

The expansion coefficients Zr(/5) may be computed by inverting eq. ( |2.9| ) using eq. 
1 



dvfii-ra) {v,w) exp [N f3 {v ■ W* + V* ■ w)] 



(A.5) 



d{l-m) ZoiP) 

The Z(/;m) iP) are actually independent of w and are of order I + m, i.e. = 0{P^~^"^) 

The factor Zq{I3) normalizes 2;(o;o) to 1 and can be explicitly computed 



where /n(x) is the Bessel function In{x) = J^^d9/{27i) exp[m6' + xcos6']. 
The Z(i.rn){P) can also be explicitly computed and read 



(A.6) 
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/iv(2/3iV) 



N . 2N^ . 
H 

{5N + 6) 2f3^N^ {7N + 12) 



+ 



(A^ + 1)^ {N + 2){N + 3) (A^ + 1)^ (A^ + 2) (A^ + 3) (A^ + 4) 
2/3" A^s (21A^3 ^ ^^3^2 ^ 214A^ + 120) 



+ . . . 



(A^ + iy{N + 2Y {N + 3){N + 4) (A^ + 5) 

N+1 , , AT + l ^ 
^(^^y = - W + + ^) "(^^^) • 

Here, we have displayed -2(i;o) to order eleven. From this expansion, one can construct 
the ordinary /3 expansions of the free energy and mass gap. 



After the expansion of each link as in eq. (|2.9| ), one has to perform the integral 
over the vectors and the link angles. For the vectors it is convenient to make use of 
the relations 



/(1;0) /(1;0) - y ^^Y^Tiy •^(2;0) , 
/(1;0)/(0;1) = fil;l) + fm , 

/(i;0)/(i;i) = ^\/(^r + l)(N + ^)^^'-''^ + ' ^^-^ 



/(o;i)/(2;i) - + 2){N + 3) ^^'-^'^ + (AT + l)(Ar + 2) ^^^^^^ " 

Appendix B: Results for the Free Energy and the Mass Gap 

We have computed the free energy to order 10 and the mass gap to order 4 in an 
^-exact strong-coupling character expansion. The result for the free energy JF is 

2 In ZM + hI =in ^1 + d„„,^ - ^) 
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q2 

^ / 2Ar 2 ^2 ^ 



e ^ ^ + 47r J 



+^(2;o)' ^ — 7- + z^-r- ^(2;o)^ 



10 / ^ 2 , / \ 8 

^1 NTi"^^'-''^ "^^'-''^ [{2n-9) {9 -in) ~ {2n + 9) {9 + in) ) '^'-''^ '^'-''^ " 
For the mass gap up to order 4 in we obtained the result 

a • 171(9, (3) — 

N ^ ( na Kb \ 
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Figure Captions 

Figure 1. A Lattice Link Configuration which Produces the Most General Plaquette- 
Variable Configuration for Open Boundary Conditions. 

Figure 2a. The Free Energy Per Unit Volume ^ as a Function of 6. Shown here 
is the order-by-order case of CP^ with P ~ 0.6. 

Figure 2b. The Free Energy Per Unit Volume ^ as a Function of 6. Shown here 
is the case of CP^ with /3 = 0.7. 

Figure 3. The Phase Diagram for the CP^ Model as Determined by the Tenth 
Order Free Energy Result. 
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Figure 4. The Fourth-Order Result for the Mass Gap am as a Function of (3 for 
Various Values of 9. 
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